Analytic Solution of the Pion-Laser Model 
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Brooding over bosons, wave packets and Bose - Einstein correlations, we find that a generaliza- 
tion of the pion-laser model for the case of overlapping wave-packets is analytically solvable with 
complete n-particle symmetrization. Multi-boson correlations generate momentum- dependent ra- 
dius and intercept parameters even for static sources. Explicit multiplicity dependence of exclusive 
correlations and spectra is found. The HBT radii are reduced for low values and enlargened for high 
values of the mean momentum in the rare gas limiting case. 
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Introduction. The study of the statistical properties of quantum systems has a long history with important recent 
developments. In high energy physics, quantum statistical correlations are studied in order to infer the space-time 
dimensions of the elementary particle reactions. In high energy heavy ion collisions hundreds of bosons are created 
in the present CERN SPS reactions when Pb + Pb reactions are measured at 160 AGeV laboratory bombarding 
energy. At the RHIC accelerator, to be completed by 1999, thousands of pions could be produced in a unit rapidity 
interval ^,^] . If the number of pions in a unit value of phase-space is large enough these bosons may condense into the 
same quantum state and a pion laser could be created Similarly to this process, when a large number of bosonic 
atoms are collected in a magnetic trap and cooled down to increase their density in phase-space, the bosonic nature 
of the atoms reveals itself in the formation of a Bose-Einstein condensate Q], a macroscopic quantum state. Such a 
condensation mechanism may provide the key to the formation of atomic lasers in condensed matter physics and to 
the formation of pion lasers in high energy particle and heavy ion physics, reviewed recently in refs [§-0. 

The density matrix of a generic quantum mechanical system is 



E 

n=0 



PnPn, (1) 



Here the index n characterizes sub-systems with particle number fixed to n, the multiplicity distribution is prescribed 
by the set of {pn}^L ' normalized as Y^=oPn = 1- The density matrixes are normalized as Trp = 1 and Tr p n = 1, 
where 

p n = Jdai...da n p n (ai, a n ) \ai, a n ){a\, a n \ (2) 

and the states |ai, a n ) denote properly normalized n-particle wave-packet boson states. 
A wave packet creation operator is 

t f d 3 p -<B=2£.-it i (p-« i ) +iu ( P xt-u) 
J {Tra z )4 

where on = 7rj, o^, tj) refers to the parameters of the wave packet i: the center in space, in momentum space, the 
width in momentum space and the production time, respectively. For simplicity we assume that all the wave packets 
are emitted at the same instant and with the same width, on — (7^, £j, a, to). 
The n boson states, normalized to unity, are given as 

I a u ... , a n ) = (^flitXiM) a& ... (4) 
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Here denotes the set of all the permutations of the indexes {1, 2, ...,n} and the subscript sized CTj denotes the 
index that replaces the index j in a given permutation from er^™- 1 . 

Solution for a New Type of Density Matrix There is one special density matrix, for which one can overcome the 
difficulty, related to the non- vanishing overlap of many hundreds of wave-packets, even in an explicit analytical manner. 
Namely, if one assumes, that we have a system, in which the emission probability of a boson is increased if there is 
an other emission in the vicinity: 



1 " / - \ 

V ' 1=1 \ CT (n) k=l / 



p n (a 1 ,...,a n )= I I />,(",){ > I I [a,,\o, T , ) j (•'>) 

The coefficient of proportionality, A/"(n) , can be determined from the normalization condition. The density matrix of 
eq. (|B|) describes a quantum-mechanical wave-packet system with induced emission, and the amount of the induced 
emission is controlled by the overlap of the n wave-packets ||, yielding a weight in the range of [1, n\]. Although it is 
very difficult numerically to operate with such a wildly fluctuating weight, we were able to reduce the problem j^] to 
an already discovered "ring" - algebra of permanents for plane- wave outgoing states . 

For the sake of simplicity we assume a non-relativistic, non-expanding static source at rest in the frame where the 
calculations are performed. 

Pi(a) = p x (0 p p (n) S(t - t ), 

= m L» exp(-e 2 /(2i? 2 )), 

ft»W - 7Z ^,3 exp(-^ 2 /(2mT)), (6) 

and a Poisson multiplicity distribution for the case when the Bose-Einstein effects are negligible: 

p°) = Hiexp(-no), (7) 
n\ 

This corresponds to the very rare gas limit, and completes the specification of the model. The plane- wave model, to 
which the multi-particle wave-packet model was reduced in ref. [|| , can be further simplified [|| to a set of recurrence 
relations with the help of the so-called "ring-algebra" discovered first by S. Pratt in ref. ||. The probability of finding 
events with multiplicity n, as well as the single-particle and the two-particle momentum distribution in such events 
is given as 

Pn = U) n I ^ Vk) , (8) 



\k=0 



N M {kl) = J2—G t (l,l), (9) 



n 1—1 

7Vi n) (kx,k 2 ) = VV [G m (l, l)Gf,_ m (2 > 2)+ 

— — 

1=2 m=l n 

+G m (l,2)G i _ m (2,l)], (10) 

where u> n = p n /po and 

G n (i,j) = Uq /i„exp(-a„(k 2 +k 2 ) +g„k i k :) ). (11) 
Averaging over the multiplicity distribution p n yields the inclusive spectra as 
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G(l,2) = 2), 



(12) 



N 1 (k 1 ) = J2p n N[ n \k 1 ) = G(l,l), 

n=l 

2V 2 (ki,k 2 ) = G(l, 1)G(2, 2) + G(l, 2)G(2, 1). 
An auxiliary quantity is introduced as 



G„ = - /d J k 1 G n (l,l) = h n - 



'o 



n \ 2a r , 



On 



With the help of the notation 



4 = cr 2 + 2mT, 



mT 

2 2 

(7 (Jrp 



(13) 
(14) 

(15) 
(16) 



the recurrence relations that correspond to the solution of the ring-algebra HU are obtained (|] for the case of the 
multi-particle wave-packet model. These correspond to the pion laser model of S. Pratt when a replacement R — > R e 
and T — ► T e = a\j[2rn) is performed. 

Let us introduce the following auxiliary quantities: 



7± = -5 (1 + x ± >/T+2s) 



r>2 2 
K e (T T 



(17) 



The general analytical solution of the model is given through the generating function of the multiplicity distribution 



where G„ is 



G(z)=£pn* r 
n=0 



G n — 



exp ^G„(z"-1) , 



7+ - 7- 



together with the general analytic solution for the functions G ra (l,2): 



G„(l,2)=j n e 



7+ ki-75k 



b n = 



1 7+ - 7- 
4 7? - 7- 



(18) 
(19) 

(20) 
(21) 



The detailed proof that the analytic solution to the multi-particle wave-packet model is indeed given by the above 
equations is described in ref. ||. 

The representation of eq. (|l^) indicates that the quantities G„-s are the so called combinants jll] [l3| of the 
probability distribution of p n and in our case their explicit form is known for any set of model parameters, as given by 



eqs. (19 17|Jl6| ) . The resulting multiplicity generating function does not correspond to discrete probability generating 
functions in ref flij , we have found a new type of probability generating functions. 

One can prove fl8J, that the mean multiplicity (n) = Y^=i n Pn = Y^aLi The large n behavior of nC n depends 

3 n 

on the ratio of no/7+, since for large values of n, we always have (7-/7+)^ << 1. The critical value of hq is 



1 + x + Vl + 2,t 



(22) 
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If uq < n c , one finds lirrin^^nCn = and (n) < oo, if uq > n c one obtains Unrin^oonCn — oo and (n) = oo, 
finally, if uq = n c one finds Zim n _,oo^Cn — 1 an d (n) = oo. The divergence of the mean multiplicity (n) is related to 
condensation of the wave-packets to the wave-packet state with n — 0, i.e. to the wave-packet state with zero mean 
momentum ||, if no > n c . 

The multiplicity distribution of eq. (|l8| ) is studied at greater length in ref. [0. 

Rare gas limiting case. Large source sizes or large effective temperatures correspond to the x » 1 limiting case, 
where the general analytical solution of the model, presented above, becomes particularly simple and the exclusive 
and inclusive spectra and correlation functions can be obtained analytically to leading order in 1/x << 1. From 
eq. ([ll|) one obtains that 



G„(l,2) = j„exp 



'la- 



-04 



*4) 



3n 



2n 

2 

x 



Ak 2 , 

(„-i) 



(23) 
(24) 



where Ak = ki — k2. We can see from eq. (p3|), that the higher order corrections will contribute to the observables 
with reduced effective temperatures and reduced effective radii. Eq. (|24|) indicates, that the leading order result 
for the combinants in the x » 1 limiting case is C\ with the first sub-leading correction given by C%. Thus, the 
probability distribution can be considered in the rare gas limiting case as a Poisson distribution of particle singlets 
with a sub-leading correction that yields a convolution of Poisson-distributed doublets. 

The very rare gas limiting case corresponds to keeping only the leading n = 1 order terms in the above equations. 
The multiplicity distribution is a Poisson distribution with (n) = no and no influence from stimulated emission. The 
momentum distribution is a static Boltzmann distribution, and the exclusive and inclusive momentum distributions 
coincide ^ . The leading order two-particle Bose-Einstein correlation function is a static Gaussian correlation function 
with a constant intercept parameter of A = 1 and with a momentum - independent radius parameter of i?* — R e [§ ■ 

This wave-packet model exhibits a lasing behavior in the very dense Bose-gas limit, which corresponds to an 
optically coherent behavior, characterized by a vanishing enhancement of the two-particle intensity correlations at 
low momentum, C(ki, k2) = 1, a case which is described in greater details in ref. B. 

The probability generating function yields the following leading order multiplicity distribution: 



Pn 



— - exp(-n ) 
n! 



1 



n(n — 1) — n^ 
2(2x)2 



(25) 



The mean multiplicity, the factorial cumulant moments of the multiplicity distribution, the inclusive and exclusive 
momentum distributions were evaluated by keeping only the leading order terms in 1 jx in ref. || . The two-particle 
exclusive correlation functions can also be evaluated by applying a Gaussian approximation to the leading order 
corrections in the x » 1 limiting case: 



>(]£!, k 3 ) 



JVj n) (ki,k 3 ) 



n(n - 1) 7v (n) 



(ki)JVi (n) (k 2 ) 



=1 + A K exp(-i? 2 C s Ak 2 - R^ Ak 2 



(26) 



where K = 0.5(ki +k 2 ), the side and outwards directions are introduced utilizing the spherical symmetry of the source 
as Ak s = Ak - K(Ak ■ K)/(K ■ K) and Ak Q = K(Ak ■ K)/(K • K), similarly to refs. |L5|]. The mean-momentum 
dependent intercept and radius parameters are 

-2\ 1 



A K =1 



(2x)i 



l-2( 5 / 2 >expf-^; 



(27) 
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Rk.s-Rl + — — 

(2x) 



(n + 2)R 2 e 



R\ - V2exp 
2 



K 2 



2 
(Tj- 



(28) 



n K 2 / K 2 



Rk,o= R k,s + ~r —4 ex P ( _ ~~T ) (29) 

Thus the symmetrization results in a momentum - dependent intercept parameter Ak that starts from a Ak=o < 1 
value at low momentum and increases with increasing momentum. Already in the first paper about the pion laser 
model, ref. ||, a reduction of the exact intercept parameter was observed and interpreted as the onset of a coherent 
behavior in the low momentum modes. First, a partially coherent system is created, characterized by Ak < 1, and if 
the density of pions is further increased, one finds a fully developed pion-laser with Ak = 0, see ref. Q for analytic 
considerations. 

Observe that the radius parameter at low mean momentum decreases while the radius parameter at high mean 
momentum increases, as compared to R e . The radius parameter of the exclusive correlation function thus becomes 
momentum- dependent even for static sources! This effect is more pronounced for higher values of the fixed multiplicity 
n, in contrast to the momentum dependence of Ak that is independent of n. 

Last, but not least, a specific term appears in the two-particle exclusive correlation function, that contributes only 
to the out direction, which, in case of spherically symmetric sources, may be identified with the direction of the mean 
momentum |l5| . This directional dependence is related only to the direction of the relative momentum as compared 
to the direction of the mean momentum, and does not violate the assumed spherical symmetry of the boson source. 



The effect vanishes both at very low or at very high values of the mean momentum K, according to eq. (29). The 
top, middle and bottom panel of Figure [l] indicates the momentum dependent Ak intercept parameter, the Rk,s 
and i?K,o radius parameters for a fireball with R = 11 fm, T = 120 MeV. The pions are assumed to be described 
by wave-packets with spatial widths of o x = 2 fm, and events with fixed identical pion multiplicity of = 600 are 
selected for the evaluation of the correlation function. For this set of parameters, the enhancement of Ak, Reside and 
Rvi,out is hardly noticeable at high momentum, but their small decrease at low momentum is clear. One may consider 
a small cold pionic system with a few large wave-packets only, by inserting R — 4 fm, T = 10 MeV, = 3 and 
o x — 5 fm to eqs. (|27|-p9|). This source could correspond to heavy ion collisions in the 30 MeV A energy domain |l5| ] 
characterized by an effective radius R e = 4.5 fm and effective temperature of T e = 15.6 MeV jus). In this case, the 
directional dependence of the radii and the enhancement of the radius parameters at high momentum as compared 
to R e becomes significant not only analytically but numerically as well. 

Highlights : In this Letter a consequent quantum mechanical description of multi-boson systems is presented, using 
properly normalized projector operators for overlapping multi-particle wave-packet states describing stimulated emis- 
sion of bosons. Our new analytic result is that multi-boson correlations generate momentum- dependent radius and 
intercept parameters even for static sources, as well as induce a special directional dependence of the correlation func- 
tion. The effective radius parameter of the two-particle correlation function is reduced for low values and enlargened 
for large values of the mean momentum in the rare gas limiting case, as compared to the case when multi-particle 
symmetrization effects are neglected. For extended, hot and rare gas of a few hundred pions, the reduction of the 
radius parameters at low momentum is found to be the most apparent effect. The directional dependence of the radius 
parameters and the enhancement of the radii at high momentum is characteristic for a small, cold pion gas with only a 
handful of particles in it. These results can be understood qualitatively by an enhancement of the wave-packets in the 
low momentum modes, due to multi-particle Bose-Einstein symmetrization effects, as the system starts to approach 
the formation of a laser, characterized by the appearance of partial optical coherence in the low momentum modes. 
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Our results explicitly depend on the multiplicity, providing a theoretical basis for event-by- event analysis of heavy 
ion data. 
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FIG. 1. Multi-particle symmetrization results at low K in a momentum-dependent reduction of the intercept parameter Ak, 
the side- wards and the outwards radius parameters, -Rk.s and Rk,o from their static values of 1 and R e , respectively. The 
enhancement of these parameters at high momentum is hardly noticeable for large and hot systems. 
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